Nanoscale magnetic field mapping with a single spin scanning probe magnetometer 
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We demonstrate quantitative magnetic field mapping with nanoscale resolution, by applying a 
lock-in technique on the electron spin resonance frequency of a single nitrogen- vacancy defect placed 
at the apex of an atomic force microscope tip. In addition, we report an all-optical magnetic imaging 
technique which is sensitive to large off- axis magnetic fields, thus extending the operation range of 
diamond-based magnetometry. Both techniques are illustrated by using a magnetic hard disk as a 
test sample. Owing to the non-perturbing and quantitative nature of the magnetic probe, this work 
should open up numerous perspectives in nanomagnetism and spintronics. 



The ability to map magnetic field distributions with 
high sensitivity and nanoscale resolution is of crucial im- 
portance for fundamental studies ranging from material 
science to biology, as well as for the development of new 
applications in spintronics and quantum technology [l]- 
[3]. In that context, an ideal scanning probe magnetome- 
ter should provide quantitative magnetic field mapping 
at the nanoscale under ambient conditions. In addition, 
the magnetic sensor should not introduce a significant 
magnetic perturbation of the probed sample. 

Over the last decades, different roads have been taken 
towards ultra-sensitive detection of magnetic fields in- 
cluding superconducting quantum interference devices 
(SQUIDs) [1], semiconductor-based Hall probes [1] and 
optical magnetometers [2^. Even though extremely high 
sensitivity has been achieved with these devices, their 
spatial resolution remains limited at the micron-scale. 
Prominent approaches to reach nanoscale resolution are 
scanning-tunneling microscopy [4l, mechanical detection 
of magnetic resonance [5 , nanoSQUIDs [6 , X-ray mi- 
croscopy [7 and magnetic force microscopy (MFM) [HI. 
Since the latter technique operates under ambient con- 
ditions without any specific sample preparation, it is 
now routinely used for mapping magnetic field gradients 
around magnetic nanostructures. However, besides in- 
troducing an inevitable perturbation of the studied mag- 
netic sample owing to the intrinsic magnetic nature of 
the probe |7], MFM does not provide quantitative infor- 
mation about the magnetic field distribution. 

Here we follow a recently proposed approach to mag- 
netic sensing based on optically detected electron spin 
resonance (ESR) [TOl. It was shown that this method ap- 
plied to a single nitrogen- vacancy (NV) defect in diamond 
could provide an unprecedented combination of spatial 
resolution and magnetic sensitivity under ambient con- 
ditions [6l [TTJ |T2] [14] . The principle of the measurement 



is similar to the one used in optical magnetometers based 
on the precession of spin-polarized atomic gases [2 . The 
applied magnetic field is evaluated by measuring the Zee- 
man shifts of the NV defect spin sublevels. In this article 
we demonstrate quantitative magnetic field mapping with 
nanoscale resolution, by applying a lock-in technique on 
the ESR frequency of a single NV defect placed at the 
apex of an atomic force microscope (AFM) tip. In addi- 
tion, we report an all-optical magnetic imaging technique 
which is sensitive to large off- axis magnetic fields, and re- 
lies on magnetic-field-dependent photoluminescence (PL) 
of the NV defect sensor, induced by spin level mixing. 

As depicted in Fig. [TJa), the scanning probe magne- 
tometer combines an optical confocal microscope placed 
on top of a customized tuning-fork-based AFM (At- 
tocube Systems), all operating under ambient conditions. 
The atomic-sized magnetic sensor was first engineered by 
grafting a diamond nanocrystal hosting a single NV de- 
fect at the end of the AFM tip HH [E]. After func- 
tionalization, a PL raster scan of the AFM tip shows an 
isolated spot at its apex (Fig.jljb)), which corresponds to 
the collected light emitted by a single NV defect, as ver- 
ified using photon correlation measurements (Fig. TJc)). 

The magnetic response of the probe was first charac- 
terized by sweeping the frequency of a driving microwave 
(MW) field through the spin resonance while recording 
the NV defect PL. Owing to spin dependent PL, ESR 
appears as a drop of the PL signal (Fig. [TJd)) [18]. The 
spin Hamiltonian of the NV defect reads as [6] 



n = hDSl + hE{Sl - Si) + h-feB ■ S 



(1) 
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where h is the Planck constant, D and E are the zero- 
field splitting parameters, z is the NV defect quantization 
axis, 7e the electron gyromagnetic ratio, and B the local 
magnetic field applied to the NV defect electron spin S = 
1. For magnetic fields smaller than roughly 10 mT, i.e. 
such that hjeB <C hD^ the quantization axis is fixed by 
the NV defect axis i tself and the ESR frequencies are 
given hy ur = D ± Y^^lTe^^^T^Trp^"^^^ where Bz is the 
magnetic field amplitude along the NV axis. The spin 
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FIG. 1: (a) -Simplified scheme of tfie scanning probe magne- 
tometer. A 20 nm diamond nanocrystal fiosting a single NV 
defect is grafted at the end of an AFM tip. A confocal mi- 
croscope placed on top of the tip allows us both to excite and 
collect the NV defect spin-dependent PL. A microwave field 
(MW) is generated by an antenna approached in the vicin- 
ity of the NV defect. (b)-PL raster scan of the AFM tip 
showing a bright PL spot at its apex. White dashed lines 
indicate the shadow of the tip. (c)- Second-order autocorre- 
lation function g^"^^ (r) recorded from the bright emission spot 
using a Hanbury Brown and Twiss interferometer. A strong 
anticorrelation effect, ^^^'*(0) ^ 0.1, is observed at zero de- 
lay. More than 1 x 10^ photons per second can be detected 
from the single NV defect when excited at saturation, with a 
signal-to-background ratio greater than 10. (d)-Optically de- 
tected ESR spectra recorded for different magnetic field mag- 
nitudes. For these experiments, the ESR contrast is C ^ 12%, 
the linewidth Aesr ~ 9 MHz, and the average rate of de- 
tected photons 71 ^ 6 X 10^ cnt/s, with an optical pumping 
power Vopt = 300 /iW. The zero- field splitting parameters are 
D = 2.87 GHz and ^ = 5 MHz. 



orientation of the magnetic sensor was determined by 
recording ESR frequencies as a function of the orientation 
and the magnitude of a calibrated magnetic field [15J. 

The shot-noise-limited sensitivity to d.c. magnetic field 
?7b is linked to the minimum detectable magnetic field 
along the NV axis B^^^ during an acquisition time At 
through the relation [19j 

r]B{T/VRi) = 5™VAt ^ - X , (2) 



where 7Z is the rate of detected photons, C the ESR con- 
trast and Aesr the associated linewidth (FWHM). From 
the optically detected ESR spectra shown in Fig. [ijd), 
we infer r]B ^ ^0 /iT/V^Hz. We note that this sensitiv- 
ity could be significantly enhanced by using pulsed-ESR 
techniques [HI [19] or by engineering a single NV defect 
in an ultra-pure diamond nanostructure [20 . 

The performances of the NV scanning probe magne- 
tometer were characterized by mapping the magnetic 
field distribution created by a commercial magnetic hard 
disk. The magnetic bit geometry of this simple test sam- 
ple is depicted in Fig. [2|a). As a first experiment, the 
magnetic sensor was approached at roughly 250 nm from 
the sample surface. At such a distance, the random orien- 
tation of the magnetic bits provides non-trivial magnetic 
field patterns with a typical magnetic field magnitude 
in the mT range [15]. The simplest way to map such 
a magnetic field distribution is borrowed from magnetic 
resonance imaging [6 . A MW field was applied with a 
fixed frequency corresponding to the zero-field ESR fre- 
quency. Iso-magnetic field images were then recorded by 
monitoring the NV defect PL while scanning the hard 
disk |6 . Dark areas corresponding to zero magnetic field 
lines can be observed (Fig.|2|b)). However, the image ex- 
hibits a varying contrast related to inhomogeneous back- 
ground luminescence from the sample. This limitation 
can be overcome by measuring the difference of PL for 
two fixed MW frequencies applied consecutively at each 
point of the scan (Fig. |2|c)). In this dual iso-magnetic 
field image, the contrast is significantly improved since 
background luminescence from the sample is suppressed. 

Even with two iso-magnetic field lines recorded at 
once, the information remains incomplete. To determine 
the full magnetic field distribution, an ESR spectrum 
could be measured for each pixel of the scan. However, 
since a few seconds are required to record such a spec- 
trum with a reasonable signal-to-noise ratio, this method 
would be extremely slow. To circumvent this limitation, 
a lock-in method was developed in order to track the shift 
of the ESR frequency while scanning the magnetic sam- 
ple [21 . This technique allows us to measure the strength 
of the magnetic field along the NV axis over the full scan- 
ning area with a typical sampling time of 110 ms [15 . As 
shown in Fig. |2|d) , the complete magnetic field map is 
in clear agreement with iso-magnetic field experiments. 
In addition, we note that the recorded pattern can be 
qualitatively reproduced by numerical simulations, both 
in shape, size, and magnetic field amplitude [15]. This 
confirms the reliability of NV-based scanning probe mag- 
netometry for imaging non trivial magnetic field distri- 
butions at the nanoscale within a reasonable data acqui- 
sition time. 

From the pixel-to-pixel noise in Fig. |2|d), an uncer- 
tainty SB ^ 27 /iT is deduced, which corresponds to a 
sensitivity r^^^ ^ 9 /nT/y/Rz, given the 110 ms sampling 
time [15] . This value is in good agreement with the shot- 
noise-limited value calculated using Eq. ([2|. Regarding 
spatial resolution, the magnetic field is experienced by 
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FIG. 2: Measurements of the magnetic field distribution created by a commercial magnetic hard disk for a probe-to-sample 
distance d — 250 nm. All images correspond to 150 x 138 pixels, with a 13-nm pixel size. (a)-Geometry of the magnetic 
bits. (b)-Single iso-magnetic field image. Dark areas correspond to a null \Bz\ component and the red arrow indicates the NV 
defect axis, along which the magnetic field projection is measured. (c)-Dual iso-magnetic field image recorded by measuring 
the PL difference for two fixed MW frequencies applied consecutively. Dark (resp. bright) areas correspond to \Bz\ = 
(resp. \Bz\ = 0.5 mT). (d)-Complete magnetic field distribution of the same area recorded by using the lock-in method. The 
panel on the right shows a line-cut taken along the white dashed line. For images (b), (c) and (d), the acquisition time per 
pixel is 5 ms, 60 ms and 110 ms respectively. 



the NV defect electron spin wavefunction, resulting in a 
subnanometric probe volume. Experimentally the spa- 
tial resolution is rather limited by the pixel size (13 nm) 
and the positioning accuracy. The 5-map depicted in 
Fig. [2|d) is thus assumed to closely follow the actual 
field distribution - within the noise SB - with an overall 
spatial resolution of a few tens of nm. The line-cut in 
Fig. |2|d) further illustrates the ability to follow abrupt 
variations, up to 0.4 mT 5B) between two adjacent 
pixels, corresponding to a field gradient of 3 x 10^ T/m. 
We emphasize the fact that no other instrument is able 
to quantitatively map magnetic field distributions with 
such a resolution. 

Yet, many material science studies are interested in 
the sample magnetization M rather than the stray field 
B. In this context, the spatial resolution is limited not 
only by the probe volume but also by its distance to 
the sample. Here for instance, simulations indicate that 
the magnetic bits could be resolved by using a probe- 
to-sample distance smaller than 40 nm (Fig. [sja)) (15j . 
However, the magnetic field amplitude then reaches sev- 
eral tens of mT, including a component B± orthogonal 
to the NV defect axis [15 . In such a situation, the quan- 
tization axis becomes determined by the local magnetic 
field rather than the NV defect axis. The eigenstates 
of the spin Hamiltonian are therefore given by super- 
positions of the rris = and mg = ±1 spin sublevels, 
both in the ground and excited states [22 . As a result, 
optically-induced spin polarization and spin-dependent 
PL of the NV defect become inefficient, and the contrast 
of optically detected ESR vanishes. Magnetic field sens- 
ing using the previous techniques is thus not accessible in 
this high off-axis magnetic field regime (> 10 mT). How- 
ever, it was also shown that the NV defect PL drops by 
roughly 30% when the magnitude of B± increases from 
to approximately 30 mT [22] [23j . This property can be 



exploited to perform a direct, microwave-free, magnetic 
field imaging using the single NV defect sensor. 

Figure [sj^b) shows a PL image recorded while scanning 
the magnetic sample without applying any microwave 
field. In this experiment, an AFM feedback loop is ap- 
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FIG. 3: (a)-Simulation of the off-axis magnetic field compo- 
nent B± for a random distribution of bit magnetizations and a 
probe-to-sample distance d = 30 nm. (b)-PL image recorded 
with the NV scanning probe magnetometer operating in tap- 
ping mode, without applying any microwave field. The image 
corresponds to 200 x 200 pixels, with a 8-nm pixel size and a 
20-ms acquisition time per pixel. (c)-Line-cut corresponding 
to the white dashed line in (b). A possible configuration of 
the bit magnetization is schematically depicted in (d). 
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plied to maintain the mean probe-to-sample distance to a 
constant value. Comparison with simulations (Fig.jsj^a)) 
suggests that this distance is around 30 nm, limited by 
the diamond size and its position on the tip as well as 
by the tip oscillation. The PL image exhibits dark ar- 
eas with a contrast greater than 20% which reveal the 
tracks and the bits of the magnetic hard disk. We note 
that the contrast depends on the field strength, which 
itself depends on the number of bits that are consecu- 
tively parallel. For instance, from the known bit-to-bit 
spacing of 55 nm, we deduce that the line-cut in Fig.jsj^b) 
corresponds to a sequence 001100, where and 1 denote 
the two possible bit magnetizations (Fig.[3]^b)-(c)). Each 
bit inversion gives rise to a stray field that appears dark 
on the PL image. We thus demonstrate a mode of NV- 
based magnetic field sensing that does not require the 
use of microwave and is sensitive to the off- axis magnetic 
field magnitude. 

Summarizing, quantitative nanoscale magnetic field 



mapping with a sensitivity below 10 /iT/ Vllz was demon- 
strated using NV-based scanning probe magnetometry. 
Moreover, we reported a microwave- free magnetic imag- 
ing technique that extends the range of application of 
NV-based magnetometry to large off-axis magnetic fields. 
Improvement of the sensitivity using pulsed-ESR tech- 
niques [19 , together with a better control of the NV-to- 
sample distance, should enable imaging of few-spins sys- 
tems, which would be of great interest for fundamental 
studies in nanomagnetism and spintronics. 
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I. SUPPLEMENTARY METHODS 



A. Experimental setup 

The experimental setup combines an atomic force microscope (AFM) and a confocal optical microscope (AFM/CFM, 
Attocube System) as shown in Fig. |4|. Two sets of piezoelectric actuators and scanners are used for positioning the 
sample and the AFM tip. The AFM cantilever is an Akiyama-Probe (Nanosensors) corresponding to a quartz tuning 
fork equipped with a micromachined silicon cantilever. The cantilever ends with a sharp tip whose nominal curvature 
radius is 15 nm. A nanodiamond hosting a single nitrogen- vacancy (NV) defect is placed at the end of the tip as 
explained in the next section. 

A continuous laser source operating at 532 nm wavelength (Spectra-Physics Excelsior) is tightly focused at the end 
of the AFM tip through a xlOO microscope objective with a 0.9-numerical aperture and a 1-mm working distance 
(Olympus, MPLFLNIOOX). The photoluminescence (PL) emitted by a single NV defect placed at the apex of the tip 
is collected by the same objective and spectrally filtered from the remaining pump light using a dichroic beam splitter 
and a band-pass filter (Semrock, 697/75 BP) . Following standard confocal detection scheme, the collected light is 
then focused onto a 50-/im diameter pinhole and directed to silicon avalanche photodiodes (Perkin-Elmer, SPCM- 
AQR-14) operating in the single-photon counting regime. Independently of the AFM scanners, precise positioning 
of the NV defect at the focus of the confocal microscope is achieved by scanning the laser beam with a fast steering 
mirror (Newport, FSM-300-02) combined with a pair of telecentric lenses. The focal length is adjusted by using the 
piezoelectric positioner of the AFM tip. 

The unicity of the emitter was checked by using a Hanbury Brown and Twiss (HBT) interferometer consisting 
of two avalanche photodiodes placed on the output ports of a 50/50 beamsplitter. The HBT detection system is 
associated with a fast multichannel analyser (Picoquant, TimeHarp 300) for recording the histogram of time delays 
between two consecutive single-photon detections. After normalization to a Poissonnian statistics, the recorded 
histogram is equivalent to a measurement of the second-order autocorrelation function g^'^\r) |T|. The observation 
of an anticorrelation effect at zero time delay, ^^^^(0) < 0.5, is the signature that a single NV defect is addressed (see 
Fig. 1(c) of the main text). 

For the electron spin resonance (ESR) measurements, a microwave synthesizer (Rohde & Schwarz, SMR20) is 
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FIG. 4: Schematic of the experimental setup, showing both the confocal optical microscope (upper part) and the tuning-fork 
based AFM microscope (lower part). L stands for lens, F for filter, A/2 for half- wave plate, P for polarizer, BS for beamsplitter, 
DBS for dichroic beamsplitter, M for mirror, FSM for fast steering mirror, PH for pinhole and APD for silicon avalanche 
photodiode. Fi (resp. F2) is a 10 nm (resp. 75nm) width interference filter centered at 532 nm (resp. 697 nm). 
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followed by a power amplifier (Mini-Circuits, ZHL-42) and connected to a 20-/im-diameter copper wire directly 
spanned on the sample surface. 

Scanning probe magnetometry is performed by monitoring the spin-dependent NV defect PL while scanning the 
magnetic sample. The tip-to-sample distance is adjusted using the piezoelectric scanner of the sample. For tapping 
mode operation, usual tuning-fork based AFM operation is used with a feedback loop adjusting the tip-to-sample 
distance. 



B. Grafting a single NV defect at the end of the AFM tip 



We started from commercially available diamond nanocrystals (SYP 0.05, Van Moppes SA, Geneva). Such nanocrys- 
tals are produced by milling type-lb high-pressure high-temperature diamond crystals with a high nitrogen content 
{[N] ^ 200 ppm). The formation of NV defects was carried out using high energy (13.6 MeV) electron irradiation 
followed by annealing at 800° C under vacuum during two hours. The electron irradiation creates vacancies in the 
diamond matrix while the annealing procedure activates the migration of vacancies to intrinsic nitrogen impurities, 
leading to NV defect bonding. After annealing, the irradiated nanocrystals were oxidized in air at 550° C during two 
hours. This procedure reduces the size of the nanodiamonds [2 and leads to an efficient charge state conversion of the 
created NV defects into the negatively-charged state 0. After sonication and washing with distilled water, aqueous 
colloidal solutions of dispersed nanodiamonds, with a mean size of 20 nm, were obtained and finally deposited by 
spin-coating onto a silica cover slip. 

A nanodiamond (ND) hosting a single NV defect was grafted at the apex of the AFM tip following the method 
introduced by Cuche et at. [4, 5 . The tip was first dipped in a solution of poly-L-lysine (EMS, molecular weight 
30000-70000 u) for a few minutes. The tip as well as the sample of dispersed nanodiamonds were then mounted in 
the experimental set-up, and the NDs were optically characterized. Once a ND hosting a single negatively-charged 
NV defect was found, the sample was scanned with the AFM tip in tapping mode in a small area around the prese- 
lected ND. Electrostatic interactions between the polymer and the ND allow reproducible and robust grafting of the 
nanodiamond at the apex of the AFM tip. 



C. Real-time tracking of the ESR frequency 



The lock-in method allowing to track the shift of the ESR frequency vr while scanning the magnetic sample is 
described in Fig. [sj Two slightly different microwave frequencies Vj^ and z^_, defined hy h'± = v ^ Az//2, are applied 
consecutively and the PL difference V{h') = 5(i^+) — 5(z/_) is computed (Fig. [sja) and (b)). This signal can be 
used as an error signal for tracking the ESR frequency [21] . Indeed, if u = ur -\- S with 5 smaller than the ESR 
linewidth (FWHM), the error signal is at first order proportional to the frequency detuning (5, i.e. to the magnetic field 
variation (Fig.jsjb)). Monitoring 0(1?) thus enables quantitative magnetic field mapping by continuously tracking the 
ESR frequency and measuring the frequency detuning from point to point, with a typical sampling time of 100 ms. 
Reliability of this method is demonstrated by monitoring in real-time the magnetic field created by an AC current 
passing through a coil with a frequency up to 1 Hz, as shown in Fig. |5|c). 



D. Estimation of the magnetic field sensitivity 



Since the magnetic stray field of the investigated hard disk is a priori not known, the measured field map (Fig. 2(c) 
in the main text) cannot be directly compared to theory. However, the field is assumed to vary slowly with respect 
to the spatial sampling, i.e. the pixel size, so that high-frequency components are attributed to measurement noise. 
Therefore, we introduce a pixel-to-pixel noise to estimate the experimental error SB. If {Bij} is the matrix cor- 
responding to the magnetic field map, with 1 < i < p and 1 < j < then SB can be defined as the standard 
deviation 

= ^ " J ^ - ^ (3) 

[p-l) X q^px {q-1) 

which takes into account the pixel-to-pixel noise along both axis. For Fig. 2(c) of the main text, for which p = 138 
and q = 150, we find SB = 27 jaT. Given the acquisition time per pixel At = 110 ms, an experimental field sensitivity 
T]^^^ = SB\fKt = 9 /iT/VHz is obtained. 
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FIG. 5: Real-time tracking of the ESR frequency. (a)-Optically detected ESR spectrum S{iy) and schematic of the field- 
frequency lock-in. (b)-Error signal T>(iy) used for the feedback loop. Provided that 6 < Aesr, ^(^) = ^ x ^ at first order, 
where A is experimentally determined and corresponds to the first derivative of V{iy) evaluated at the resonance frequency ur. 
For each cycle of the tracking algorithm, T>(i') around the former ESR frequency is measured and the local magnetic field Bz 
is deduced. (c)-Real-time measurement of the magnetic field created by a 0.5 Hz AC current in a coil with a sampling rate of 
150 ms. Parameters of the feedback loop are Az/ = 5 MHz and ^ = 5 X 10^ cnts.s"^MHz"^ 



E. Determination of the NV defect axis 



The NV defect axis, characterized by the unit vector can be determined by measuring the ESR spectrum while 
applying a static magnetic field along several well-chosen directions. Denoting 6 and (j) the polar and azimuthal angles 
of the applied magnetic field in the reference frame xyz of the NV defect, the diagonalization of the electron spin 
Hamiltonian (equation (1) of the main text) yields to the relation [6] 

7D^ + 2(^1 + (2 {u^ + uj) - 5^1 - 9E^) - + ^| - ^1^2 + 9^^) 

where ui and 1^2 are the Zeeman-shifted ESR frequencies, A = Dcos20 + 2£^ cos 2(/) sin^ ^ while D and E are the 
zero- field splitting parameters of the NV ground-state electron spin. Since E <^ one can use the approximation 

For a magnetic field applied along a given direction Eq. Q gives the angle \6k\ between z and K. Repeating 
the measurement for two other directions of the magnetic field enables to completely determine z. Since the magnetic 
images presented in the main manuscript are XY maps of the magnetic field component \Bz\^ it is important to know 
z with respect to the reference frame XYZ of the laboratory. The orientation of the NV defect used in Fig. 2 and 3 of 
the main text was determined following the latter method. Fig. [6|b) shows the evolution of the measured resonance 
frequencies ui and 1^2 as a function of a static magnetic field produced by two coils in an Helmholtz configuration. 
This field was first aligned along the X-axis (filled circles) and then along the F-axis (open circles). From these 
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measurements and Eq. Q, the angles \0x\ = 37 ± 1° and |6>y | = 53 ± 1° were obtained. As |6>x| + | = 90°, one 
deduces that z hes in the XY plane. The actual orientation of the NV axis in the reference frame XY Z is illustrated 
in Fig.|6ja). The other orientation compatible with the set (|^xM^y|) was ruled out by a measurement along a third 
direction. 




CoilY Bcoil(mT) 



FIG. 6: Measurement of the NV defect orientation using Helmholtz coils. (a)-Measuring the ESR spectrum while applying a 
static magnetic field along X (resp. Y) gives access to the angle \0x\ (resp. \0y\) between the NV axis (red arrow) and the 
X-axis (resp. F-axis). (b)-ESR frequencies vi and V2 as a function of the magnetic field magnitude. The filled (resp. open) 
black circles are the experimental data for a static magnetic field applied along X (resp. Y). The angles \0x\ and \0y\ were 
estimated using Eq. (|4|. The solid lines are the theoretical curves using the full diagonalization of the spin Hamiltonian and 
the calculated values of \0x\ and considering a field oriented along X (green lines) and Y (blue lines), respectively. 



II. MODELLING OF THE MAGNETIC FIELD DISTRIBUTION ABOVE THE HARD DISK 

The magnetic sample used for the demonstration of magnetic field mapping is a piece of a commercial magnetic 
hard disk (Maxtor, 20 Gb). The geometry of the magnetic bits is depicted in Fig. [?[ AFM topography of the sample 
allows to infer the geometrical orientation of the bits since grooves imprinted by the read/write head indicate the 
direction of the tracks (Fig. [7|^a)). This information is used to determine the NV orientation with respect to the bits 
orientation. 

The hard disk was modeled as an array of uniformly magnetized parallelepipeds (the 'bits'). The dimensions of 
each 'bit' are 2a = 40 nm, 2h = 400 nm and 2c = 10 nm, along X, Y and Z, respectively, and their magnetization 
is M = iMgX depending whether the bit is at '0' or '1' (Fig. [7|b)-(c)). Mg is the average magnetization of the 
ferromagnetic material. Typical value for the studied hard disk is ijlqMs = 0.5 T where /io is the permeability of free 
space. The pitch of the array is 55 nm along X and 450 nm along Y . 





FIG. 7: Modelling the magnetic hard disk, (a)- AFM Topography of a 2 x 2 /xm area of the hard disk. (b)-Model used for the 
20 Gb hard disk. Each bit is magnetized along ±X and has dimensions 400 x 40 x 10 nm. The bits are spaced by 450 nm 
along Y (track width) and 55 nm along X. (c)- A parallelepiped uniformely magnetized along X is equivalent to a magnetic 
"capacitor" with two charged surfaces located dX X — zba. 



9 



The magnetic field distribution produced by a single bit (Fig. [7|^c)) can be calculated using magnetostatics theory. 
Outside the parallelepiped, the three components of the magnetic field read as 



Bx = 



in 



Yl ^^■'^t^^ ' 



i=±ij=±i k=±i 



{Y + jb){Z + kc) 



{X + ia) v/(xTi^)2TTrTj6)2TpTfc^ 



(5) 



By 



Stt 



In 



T-r T-r T-r ( {Z + kc) - ^{X + iaf + (F + jbf + (Z + kef 



ijk' 



(6) 



B, 



Stt 



In 



TT TT TT ( {Y+ jb) -^/{X + iaf + {Y + jbf + {Z + kef ' 

hLMikX \{Y+jb)+^w+w+(x+w+w+w, 



ijk' 



(7) 



The magnetic fields produced by each bit of the array were then summed up in order to compute the whole field 
distribution above the hard disk. Since the exact distribution of magnetic bit orientations is unknown, we compute 
the magnetic field distributions along X, Y and Z directions for a random distribution of magnetic bits (Figure [s]^ a)). 
Figure [sjb) indicates the magnetic field distribution at 250 nm from the hard disk, once projected along the axis of 
the NV defect used in the scanning-probe magnetometry experiments. The magnetic field follows rounded non trivial 
patterns with characteristics size and amplitude similar to the one observed experimentally (Fig. 2(c) of the main 
text). We note that the magnetic field distribution at another distance or along the Y axis follows a clearly different 
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FIG. 8: (a)-Calculated magnetic field distributions above a randomly written area of the hard disk. The three components 
Bx^ By and Bz (columns 1, 2 and 3, respectively) are calculated at various distances d from the surface, from 10 nm to 250 
nm (from bottom to top). The plots are 1.8 x 1.8 /xm in size, but the array used for the calculation is 9 times larger (5.4 x 5.4 
jjLVDi) in order to avoid edge effects, (b)- Simulation of the magnetic field created by a random distribution of bit magnetizations 
for d — 250 nm, once projected along the NV defect axis used in the experiments described in the main text of the article. The 
magnetic bits are oriented as shown in the AFM topography (Fig. 4(a)). 
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